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Abstract. The emergence of phase separation is investi- 
gated in the framework of a 2D t— J model by means of a 
variational product ansatz, which covers the infinite lattice 
by two types of L x L clusters. Clusters of the first type 
are completely occupied with electrons, i.e. they carry max- 
imal charge Q e — L 2 and total spin 0, and thereby form 
the antiferromagnetic background. Holes occur in the sec- 
ond type of clusters - called "hole clusters" . They carry a 
charge Qh < L 2 . The charge Qh and the number N{Qh) of 
hole clusters is fixed by minimizing the total energy at given 
hole density and spin exchange coupling a = J/t. For a not 
too small (a > 0.5) it turns out that hole clusters are occu- 
pied with an even number Qh < L 2 of electrons and carry a 
total spin 0. For increasing a the charge Qh{a) of the hole 
clusters decreases. Some points on the boundary curve can 
be extracted from Qh(a). 
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1. Introduction 

Holes play a fundamental role in our understanding of high T c superconductivity, jjj The 
parent materials like La 2 CuO^ are insulators with antiferromagnetic order and doping 
with holes (missing electrons) opens the superconducting phase. [2j 

Experimental evidence has been found in La 2 CuO^ for phase separation. |SHH El IE] 
This means that the holes in the Cu0 2 planes are not distributed uniformly but 
concentrate in "hole rich" domaines. The compound phase separates for 0.01 < 5 < 0.06 
below T ps ~ 300K into the nearly stoichiometric antiferromagnetic La^CuO^^ with 
5\ < 0.02 and Neel temperature (Tjy — 250K), and a metallic superconducting oxygen- 
rich phase La 2 CuO i+ s 2 with 5 2 ~ 0.02 and T c ~ 34K. The other evidence is related to 
the Sr doped compound La 2 _ x Sr x CuO i+1 . This compound phase separates for x < 0.03 
into the superconducting La 2 - x Sr x CuOi +11 (71 ~ 0.08) and the nonsuper conducting 
La 2 - x Sr x Cu0 4:+72 (72 ~ 0.00) phases. [7] 

Intensive studies have been performed to understand this phenomenon in models for 
strongly correlated electrons like the Hubbard - and t — J model. In the 2D t — J model 
various attempts have been made to exploit the phase diagram in the plane spanned 
by the charge density p = Q^/N and the spin coupling a = J/t (N, Qtot denote the 
total numbers of sites and electrons, respectively). 

Phase separation occurs, if both the charge density p and the spin coupling a, are 
large enough 

1 > P > Pi, a > ap(pi) • (1) 
In this regime, the ground state can be represented by a product ansatz 

4(p,N)=4 K (N e (p))l> h (N h {p,p 1 j) (2) 
of two clusters with site numbers N e (p) and Nh(p,pi), which cover the whole lattice: 

N e (p) + N h (p, Pl ) = N. (3) 

The N e (p) sites in the "electron" cluster are all occupied with electrons. The 
corresponding ground state ip e (N e (p)) is just the ground state of the 2D Heisenberg 
model with N e (p) sites. Holes occur in the "hole" cluster with N h (p, pi) sites. The 
corresponding ground state iph(N h (p, pi)) is given by the t — J model with N h (p, pi) 
sites and Qh = Nh(p, pi)pi electrons. The numbers N e (p) and Nh(p, pi) for the cluster 
sites are fixed by the total number of sites (jSJ) and the total number of electrons 

N e (p)+ Pl N h (p, Pl ) = Q tot . (4) 

From 0-dlI) one derives that the ground state energy per site e(p, a): 

e{p,a) = [e(pi,a)(l - p) + e(l,a)(p - pj] (5) 

1 - Pi 

is linear in p. This holds for i.e. in the region with phase separation. 

The focus is the boundary curve a p (pi) for phase separation. Two controversary 
points of view can be found in the literature: Emery et al.jH] suggested that in the 
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2D t — J model the phase separation curve a p (pi) starts at a p (pi = 1) = 0, which 
means that phase separation exists already for low doping (1 — pi) ~ and small values 
of a, as observed experimentally. |3] Their point of view is supported by Hellberg and 
Manousakis CHI HI C21 • 

On the other hand Putikka et al. ITl] concluded from a high temperature 
expansion, that phase separation at low doping only emerges for larger a values 
a > a p (p\ — 1) ~ 1.2, which would exclude the a regime realized in the experiments. 
This point of view is supported by DMRG calculationss on ladders (Rommer, White, 
ScalapinojT3]), variational wave functions of the Luttinger-Jastrow-Gutzwiller type 
(Valenti, Oros|f tij). Lanczos calculations (DagottofTTj llHj). 

Finally, Green's function Monte-Carlo simulations performed by Calandra et al. 
(CBS) [IH! led to a value a p (p x = 1) ~ 0.5. 

In this paper, we would like to study the mechanism of phase separation by means 
of a product ansatz with Lx L clusters of charge Q which cover the infinite lattice. Our 
method is in the spirit of the coupled cluster method (CCM) designed almost 50 years 
ago [2D] as an approximation scheme for quantum many body problems. In order to 
handle the interaction between clusters, perturbative [21] and variational [22] methods 
have been implemented. In our study of the phase diagram in the (p = Q/N,a = J/t) 
plane, we proceed as follows: We start from the ground states ip^Hpi, a ) with energies 
E^ L \Qi,a) on L x L clusters with charge Q\ = p\L 2 . These ground states can be 
computed analytically for plaquette clusters (L = 2). In Section |21 we will discuss the 
product ansatz with plaquette clusters, which minimizes the total energy at fixed charge 
density p = Q/N . Phase separation can be observed on the small plaquette cluster on 
a "microscopic scale". 

In Section 01 we extend our considerations to clusters of size Lx L. Results for the 
phase diagram obtained from a numerical calculation of ground state energies on a 4 x 4 
cluster are shown. We finally summerize and present the discussion in Section |3J 

2. Plaquette cluster in the t — J model: Phase separation "in statu 



We start from the t — J model in two dimensions which is built up as a lattice of L x L 
clusters (cf. Fig. HJfor the example of L = 2): 



with h Pi and h PuPj containing all nearest neighbour bonds (k, I) in a cluster pi or between 
neighbouring clusters (pi, pj), respectively. In general, each nearest neighbour bond (k, I) 
of the lattice either belongs to a cluster (parameters t, J = at) or connects neighbouring 
clusters (t",J" = at") and contributes with th^ k,l \a) or t"h^ k,l \a) to the respective 
intra- or inter-cluster part of the Hamiltonian ©: 



nascendi 



H = tY,h Pi {a= J/t)+t'"£h Pi!P .(a) 



(6) 



(id) 




a 
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where 

ni = Y. n ^ = H c t* c i,° ( 8 ) 

and V projects onto the subspace where the occupation numbers ni are restricted to 
zero or one. The electron spin operator is represented by at the k-th site. 

The interaction between L = 2 clusters is mediated by the two dashed links for 
each spatial dimension. It is treated in the following by a perturbation expansion in 
the hopping parameter t", which allows the hopping of electrons between neighbouring 
plaquettes. 
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Figure 1. 2D lattice structure with underlying 2x2 plaquette structure. 



To zeroth order in t", the eigenstates of the 2D t — J model decay into a product 

U\Qp) (9) 

p 

of plaquette eigenstates \Q P ) with charge Q p . The ground state is characterized by a 
minimum of the energy: 

E = 5><">(<W = E E^\Q)N{Q) (10) 

P Q=0 

where N(Q) is the number of plaquettes with charge Q = 0,1,2,3,4. E^ P \Q) is the 
ground state energy of a plaquette with charge Q. 

The numbers N(Q) are constrained by the total number of plaquettes (N/4) and 
the total charge (Qtot) : 

4 ^4 

£ N(Q) = N<*> = - , £QiV(Q)=Q tot (11) 
Q=o 4 <2=o 

The plaquette ground state energies E<p*(Q) have been computed in Appendix 
A of Ref. [23- The there obtained formulas simplify for the symmetric case we are 
considering here. 
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In the charge sectors Q = 0,1,2,4 there is no level crossing and the plaquette 
ground state energy for all a > reads 

£ (p) (0) = 5 = (12) 

E {p \l) = - 2t S = l/2 (13) 

£?W(2) = -7}(a + Va 2 + 32) S = (14) 

£(p)(4) = -3ta S = (15) 

Note that the total spin of the plaquette electrons commutes with the plaquette 
Hamiltonian and the eigenvalues of = S(S + 1) can be used to characterize the 
plaquette ground state. This plays an important role in the sectors with Q = 3 where 
the ground state changes with a: 

E { f ) {3, a) = -2t, S = 3/2 for < a < a m 



E^(3,a) = -t(a + \/x + 3 )' s=1 / 2 for a m <a<2 



E^j(3, a)= - t ( y + l) , S=l/2 for a > 2. (16) 

with a m = |(4 - v 7 !^) = 0.262... 

Note, that the Nagaoka ferromagnetic state [21] - with maximal spin S = Q/2 - 
is already visible on a 4-site plaquette with one hole, i.e. Q = 3, if a is small enough 
a < a rn = 0.262. Eigenstates with maximal plaquette spin S = Q/2 exist in all charge 
channels; the corresponding energy eigenvalues E(Q,S = Q/2) do not depend on the 
spin exchange coupling a. However, these eigenstates are not ground states except for 
the one hole case. 

Let us next look for the minimum of the total energy ([10)1 . which depends on the 
relative magnitude of the plaquette energies. 

a) In the regime, where the inequality 

£ (p) (4) + E (p) {2) < 2E (p \3) (17) 

holds, two plaquettes with charge Q = 3 are substituted by two plaquettes with 
charges Q = 2 and Q = 4, respectively. Insertion of the ground state energies [(JHJ), 
(fTBl . ([ISD ] yields the validity of (O for 

a > a(3) = -jL ~ 0.436 . (18) 



b) The inequalitiy 

E (p \0) + E {P) (A) < 2E {p) (2) (19) 

holds for 



a > a(2) = 4^2/3 ~ 3.266 (20) 
and allows the substitution of the charge Q = 2 plaquettes. 
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c) The inequality 

£ (p) (0) + E {p) {2) < 2E {p \l) (21) 

holds for 

a > 2 (22) 
and allows to substitute Q — 1 plaquettes. 

Fig. EJA) shows a schematic view of the phase diagram in the p — a plane with the 
pairs of plaquette charges, which fix the ground state energy per site according to eqs. 

(uni-tnu). 
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Figure 2. Schematic phase diagram obtained from 2x2 plaquettes (A) and 4x4 
plaquettes (B). The charges Q\,Qi that constitute the respective ground states in the 
p — a-planes and the boundaries a(2j) are given in the text. 



In each of the domains with plaquette pairs Qi, Q 2 

w c <p< w c 

where the plaquette sites iV c = 4, the ground state energy per site: 

1 -{^ (p) (Qi)(t-p)+^ ) (Q 2 )(p-f)} 



e(p,a) 



Q2-Q1 



(23) 



(24) 



is linear in p and the chemical potential p(p,a) = de/dp is constant with respect to p. 
Note that two domains (Qi, Q2) and (Q2, Q3) with a common plaquette charge Q2 have 
also a common horizontal boundary at charge density p 2 = Q%jk. At this boundary the 
chemical potential p(p, a) as function of p is discontinuous with a jump 

A(Q 2 , a) = p(p + 0, a) - p(p - 0, a) 

Q3-Q1 



Q3-Q2 Q2-Q1 (Q3-Q2XQ2-Q1 



-E®{Q 2 ) (25) 
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which vanishes for a specific value of a = a(Q 2 ). They define the vertical lines 
in the phase diagram given by (|18|h ()20|) and (|2*2*j) . This type of phase diagram 
has been introduced first by Kagan et al. [2H] for the three leg ladder in the rung 
cluster approximation. Note, that the constraints for (N^, 0{- o t) (jll|) are taken into 
account explicitely on both sides of the inequalities (fTTj) . (fTTIJ) . (|2*T]) . This procedure 
can be considered as an alternative to the usual grand canonical one, where the charge 
conservation constraint (jllj) is eliminated by a reservoir with a fixed chemical potential. 



3. Product ansatz for the ground state with large clusters 

The considerations which led to the phase diagram in Fig. EtA) based on plaquette 
clusters can easily be extended to larger L x L = N c clusters with N c sites and 
charges Q = N c , . . . , 0. We start with the ground state energies E(Q, N c , a) and their 
dependence on a. The generalization of (|24jl to a cluster with N c sites and charges 
Qi < Q2 is straightforward by replacing E^ p \Qi) — > E(Q 1 , a, N c ). Again the ground 
state energy per site is linear in p and the chemical potential is constant: 

MP, a) = - 1 - {E(Q 2 , a, N c ) - E(Q U a, N c )\ . (26) 

Two domains (Qi,Q2) {QiiQs) with a common cluster charge Q2 have a common 
boundary in the phase diagram at p 2 = Q 2 /N c . At this boundary the chemical potential 
is discontinuous with a jump 

A(Q 2 , a) = p(p 2 + 0, a) - p(p 2 - 0, a) 

= E(Q 3 )-E(Q 2 ) _ E{Q 2 )-E{Q l ) 
Q3 — Q2 Q2 — Qi 

which vanishes for a specific value of a = a(Q 2 ). The inequality 

A(Q 2 ,a)<0 for a > a{Q 2 ) (28) 

means that the clusters with charge Q 2 and ground state energy E(Q 2 ,a, N c ) can be 
substituted by clusters with charges Q 3 and Q\ and energies E(Q 3 , a, N c ), E(Qi, a, N c ). 
Here the two domains (Qi,Q 2 ) {QziQz) merge for a > a{Q 2 ) and the ground state is 
given by a cluster product ansatz with charges {Qi,Qz)- 
Let us first consider the case [case (A)] 

Qi = Qi- 1, Qs = Q2 + 1, Q2 odd, a < a{Q 2 ) . (29) 
Here the ground state is built up from clusters with charges: 

(Q 2 -1,Q 2 ) for ^^<P<^ (30) 

(Q 2 ,Q 2 + 1) for 9l <p<^. (31) 

The vanishing of the jump (J27j) in the chemical potential 

A(Q 2 , a(Q 2 )) = E{Q 2 + 1) + E(Q 2 - 1) - 2E(Q 2 ) = (32) 
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Table 1. Boundary couplings a{Q2) [from left to right cases (A), (B), (C)) where the 
jump (|27|) in the chemical potential vanishes for cases (Qi, Q2, Qz)- 
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defines the boundary a = a(Q 2 ), where the two domains (}30|) and (j3~Tf merge together: 
(Q2-I, Q2 + I) for ^_l<p<^±l anda>a(Q 2 ). (33) 

The numerical evaluation of (|32|) from the ground state energies on a 4 x 4 = 16 cluster 

- with periodic boundary conditions - yields the couplings a(Q 2 ) listed in Tabled (col. 
1,2): 

In the next step, we consider the cases [case {B)\. 

<2i = <3 2 -2, Q 3 = Q 2 + 2, Q 2 = 14,12,8,4, a < a{Q 2 ) . (34) 
Here the ground state is built up from clusters with charges 

(Q 2 -2,Q 2 ) for 9izl< p <^ (35) 

(Q 2 ,Q 2 + 2) for 9l< p <93±l. (36) 
The vanishing of the jump (}2*Tj) in the chemical potential 

A(Q 2 , a(Q 2 )) = l - [E(Q 2 + 2) + E(Q 2 - 2) - 2E(Q 2 )] = (37) 

defines again the boundary a = a(Q 2 ), where the two domains (j33j) and (|3l)j) merge 
together: 

(Q2-2, Q 2 + 2) for ^l< p <^±^ and a > a(Q 2 ). (38) 

The numerical evaluation of ()37|) from the ground state energies on a 4 x 4 = 16 cluster 

- with periodic boundary conditions - yields the couplings a{Q 2 ) listed in Table [T] (col. 

3,4): 

We are left with clusters of charge Q 2 = 10, 6, 2 [case (C)]. 

They are eliminated successively by considering the jumps (|2*7j) in the chemical 
potential for the cases in Table ^ (col. 5,6). These jumps vanish at the couplings 
a = a(Q 2 ) listed in the third column, such that the clusters with charges Q 2 can be 
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eliminated for a > ol{Q<2). The resulting phase diagram is shown in Fig. |2fB). The pairs 
of integers in the rectangular domains denote the two cluster charges which determine 
the ground state in the product ansatz. 

Let us comment the appropriate boundary conditions for the L x L cluster in the 
variational product ansatz for the ground state on the infinite lattice. A priori we are 
free in our choice of the clusters and their boundary conditions. In our opinion periodic 
boundary conditions are most appropriate for the following reasons: 

The variational product ansatz becomes exact in the limit of infinite cluster (L — > 
oo, N c = L? — > oo). In this case, the cluster energies E(Q, a, N c ) = N c e{p\ = Q/N c , a) 
are related to the ground state energies per site e{p\ = Q/N c , a) at fixed charge density. 
On finite clusters with Lx L = N c sites the ground state energies are approximated most 
accurately with periodic boundary conditions, since the interaction between the clusters 
is partly taken into account by means of the periodic boundary terms. Of course - in a 
perturbation expansion for the interaction between the clusters - the periodic boundary 
terms have to be subtracted again in each order. 



4. Summary and discussion 

In this paper we made an attempt to describe phase separation in terms of a product 
cluster ansatz for the ground state of the 2D t — J model. 

The analytic results for 2 x 2 plaquette clusters reveal already the gross features as 
depicted in the phase diagram in Fig. |21 For p and a large enough the ground state is 
built up from two types of plaquettes (Qh, Qe) = (0, 4), (2, 4). The hole clusters carry a 
charge Qh which increases with decreasing a: 



Q h = a(2) = 4W-<a 0<p<l (39) 




Q h = 2 a(3) = -|= < a < a{2) \<P<1 (40) 



where p\ = Qh/N denotes the charge density in the hole cluster. The second type of 
plaquettes is completely occupied with electrons Q e = 4. 

In the regimes ()39j) and (f4T)|) the ground state energy per site ()24j) is linear in p in 
the respective intervals. The lower bounds in the a-intervals listed in (|39|). ()40|) define 
two points on the boundary (JTJ) for phase separation: 

l~2 12 
ocp{p% = 0) = 4W-, a p (pi = -) = —=. (41) 



The numerical results - obtained from a product ansatz with 4x4 clusters lead 
to the phase diagram in Fig. El^B) similar to Fig. W[A) for 2x2 clusters. In 
each rectangular domain the cluster charges Qi,Q2 result from a minimization of the 
ground state energy per site (|24*j) at fixed charge density p (in the infinite system). 
Phase separation is observed in the upper part of Fig. EfB) with cluster charges 
(Qh,Qe) = (0, 16), (2, 16), (6, 16), (10, 16) The charge Qh of the hole clusters changes 
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with the spin exchange coupling a 



Q h = a(2) = 3.5 < a, < p < 1 (42) 

Q h = 2 a(6) = 2.2 < a < a(2), - < p < 1 (43) 

8 

3 

Q h = 6 o(10) = 2.1<a<a(6), - < p < 1 (44) 

8 

Q h = 10 a(14) = 0.6 < a < a(10), - < p < 1 (45) 



The ground state energy per site (|24jl is linear in the p intervals listed in (|42)1-(|4T 
The lower bounds in the a intervals yield 4 points on the boundary curve (0) for phase 
separation shown in Table El 



Table 2. Points on the boundary curve (JIJ for phase separation derived from a 4 x 4 
cluster. 
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The product ansatz (cf. for plaquettes) with isolated cluster ground states is 
highly degenerate. Each distribution of electron and hole clusters over the whole lattice 
leads to the same ground state energy, if we neglect the interaction between neighbouring 
clusters (cf. Fig £2). 





Figure 3. Interaction between neighbouring plaquettes with charges Q\,Q2 shown 
for a pair of nearest neighbour sites x and x' . 

Let us denote the interaction energy between neighbouring clusters with charges 
Q 2 = L 2 , Q 1 = Q h by W(Q h Q 2 ) (cf. FigEJ. If the difference 

A = W(L 2 , L 2 ) + W(Q h , Qh) - 2W(Q h , L 2 ) < (46) 

is negative, the ground state prefers phase separation in the following sense. The 
numbers N(L 2 ,L 2 ), N(Qh,Qh) of identical neighbouring clusters is maximal 
N(LM* )=2 N m N {Qh ,Q h)=2 Nm N ^ oo 

N N N N K ' 

where N(L 2 ) and N(Qh) are the numbers of electron and hole clusters, respectively. As 
a consequence the number of neighbouring clusters N(Qh, L 2 ) with different charges is 
minimal: 

^% -> for N^oo. (48) 
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Within the product ansatz (JUJ) with plaquettes of charge Q 2 = L 2 , Qi = Q h the 
interaction energies turn out to be 

W(Q 1 ,Q 2 )=J2(QuQ2\h {x ' x,) \Qi,Q 2 )= ~^QiQ 2 , (L = 2), (49) 

provided that the cluster charges are even and the total cluster spins are zero. In this 
case the spin matrix elements 

(QilSalQi) = (Q 2 \S x \Q 2 ) = (50) 

and the hopping contributions arising in (J2J) vanish for L 2 — Qh > 2. Eq. ()49|) implies 
that the inequality flfij) is valid and the ground state configuration (|4*7j) . (f4*KJ) with two 
big clusters of charge density p 2 — 1, Pi = Qh/L 2 is selected out. According to (|4*2^) -(|43 j) 
the necessary condition L 2 — Qh > 2 is satisfied for a > 0.6. 

The phase diagrams in Fig. EfA) and Fig. B) contain the whole information on 
the ground state of the system provided that a product ansatz with two L x L clusters 
and charges Q\, Q 2 is adequate. Since we know the cluster ground states in each 
charge sector from the analytical calculation in Ref. [23J for L = 2 and our numerical 
calculation for L = 4, we can also determine the hole-hole correlators: 

C ix ,y)(Qi) = (n h (0,0)n h (x,y)) - (n h (0,0))(n h (x,y)) (51) 

where 7^(0,0), Uh(x,y) count the number of holes at sites (0,0) and (x,y). 

On the 4x4 cluster with periodic boundary conditions 5 independent correlators 
can be arranged according to the distance vector (x,y). Let us look at the following 
regime, 

2.2<a<3.5, -<p<l. (52) 

8 

According to the phase diagram in Fig. EfB) the ground state contains clusters with 
charges Qi = 2, Q 2 = 16. The hole-hole correlators ([ST]) in the Q\ = 2 cluster turn out 
to be negative. However, the modulus of |Ci )0 (a)| decreases with a, whereas |C 2) o(a)|, 
|C 2) i(a)| and |C 2j 2(a)| increase with a. We interprete these increasing "long"-range 
correlations as a hint to condensation of holes for large a- values. 
Next let us look for the hole-hole correlators in the a-regime 

0.6<a<2.1, -<P<1- (53) 

8 

According to the phase diagram in Fig. EfB) the ground state contains clusters 
with charges Q± = 10, Q 2 = 16. Again all correlators are negative. The "long"- 
range correlators | C 2 ,o \ = |Ci,i|, | C*2,i | , \C 2 , 2 \ are small in comparison with the nearest 
neighbour correlator | Ci,o I for 0.6 < a < 1.5. In this regime the system prefers formation 
of hole pairs on neighbouring sites. 

On the other hand for a < 0.6 and p > 0.75, we observe rapid changes of the ground 
state with p and a as is indicated by the numerous small rectangular domains in the 
left upper part in Fig. EfB). In this regime holes are no longer confined as can be seen 
directly from the product ansatz with 2x2 plaquettes and charges Qh = 3, Q e = 4 [cf. 
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upper left part in Fig. 12(A)]. Here, the hopping matrix elements ((7j) are active already in 
first order perturbation theory and allow for the exchange of hole (Qh = 3) and electron 
(Q e = 4) clusters. In principle, the holes can now hop over the whole lattice and destroy 
thereby phase separation. The distribution of the holes can be determined only from a 
precise computation of the ground state in the low doping 5, low a regime. 

References 

[1] Bednorz J G and Miiller K A 1986 Z. Phys. 5.64 189 
[2] Anderson P W 1987 Science 235 1196 

[3] Jorgensen J D, Dabrowski B, Pei S, Hinks D G and Soderholm L 1988 Phys. Rev. B. 38 11337 
[4] Hammel P C, et. al 1990 Phys. Rev. B.42 6781 
[5] Hammel P C, et. al 1991 Physica C185-189 1095 
[6] Chou F C, et. al 1996 Phys. Rev. 5.54 572 
[7] Cho J H, et. al 1993 Phys. Rev. Lett. 70 222 

[8] Emery V J, Kivelson S A and Lin H Q 1990 Phys. Rev. Lett. 64 475 

[9] Hellberg C S and Manousakis E 1995 Phys. Rev. 5.52 4639 
[10] Hellberg C S and Manousakis E 1997 Phys. Rev. Lett. 78 4609 
[11] Hellberg C S and Manousakis E 1999 Phys. Rev. Lett. 83 132 
[12] Hellberg C S and Manousakis E 2000 Phys. Rev. 5.61 11787 
[13] Putikka W O, Luchini M U and Rice T M 1992 Phys. Rev. Lett. 68 538 
[14] Putikka W O, Glenister R L, Singh R R P and Tsunetsugu H 1994 Phys. Rev. Lett. 73 170 
[15] Rommer S, White S R and Scalapino D J 2000 Phys. Rev. B. 61 13424 
[16] Valenti R and Gros C 1992 Phys. Rev. Lett. 68 2402 

[17] Dagotto E, Moreo A, Ortolani F, Poilblanc D and Riera J 1992 Phys. Rev. B 45 10741 

[18] Dagotto E, Martins G B, Riera J, Malvezzi A L and Gazza C 1998 Phys. Rev. B 58 12063 

[19] Calandra M, Becca F and Sorella S 1998 Phys. Rev. Lett. 81 5185 

[20] Coester F 1958 Nucl. Phys. 7 421; Coester F and Kummel K H 1960 ibid. 17 477 

[21] Gros C and Valenti R 1993 Phys. Rev. 5 48 418; Senechal D, Perez D and Pioro-Ladriere M 2000 

Phys. Rev. Lett. 84 522; Senechal D, Perez D and Plouffe D 2002 Phys. Rev. B 66 075129 
[22] Potthoff M 2003 Eur. Phys. J. B 32 429; Potthoff M 2003 Phys. Rev. Lett. 91 206402 
[23] Fledderjohann F, Langari A and Mutter K H 2004 J. Phys.: Condens. Matter 16 8571 
[24] Nagaoka Y 1966 Phys. Rev. 147 392 

[25] Kagan M Y, Haas S and Rice T M 1999 Physica C 317 185 



